Effect of Coulomb interaction on the gap in monolayer and bilayer graphene 
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We study effects of a repulsive Coulomb interaction on tlie spectral gap in monolayer and bilayer 
graphene in the vicinity of the charge neutrality point by employing the functional renormalization- 
group technique. In both cases Coulomb interaction supports the gap once it is open. For monolayer 
graphene we correctly reproduce results obtained previously by several authors, e.g., an apparent 
logarithmic divergence of the Fermi velocity and the gap as well as a fixed point corresponding to a 
quantum phase transition at infinitely large Coulomb interaction. On the other hand, we show that 
the gap introduces an additional length scale at which renormalization flow of diverging quantities 
saturates. An analogous analysis is also performed for bilayer graphene with similar results. We find 
an additional fixed point in the gapless regime with linear spectrum corresponding to the vanishing 
electronic band mass. This fixed point is unstable with respect to gap fiuctuations and can not be 
reached as soon as the gap is opened. This preserves the quadratic scaling of the spectrum and 
finite electronic band mass. 

PACS numbers: 71.10.-w, 05.10.Cc 
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I. INTRODUCTION 

Monolayer (ML) and bilayer (BL) graphene are 
semimetals with an electron and a hole band. Both bands 
touch each other at two nodes. The low-energy disper- 
sion in the vicinity of these nodes is linear in ML graphene 
and quadratic in BL graphene. Exactly at the nodes both 
systems obey a chiral symmetry which reflects the sub- 
lattice symmetry of the underlying honeycomb lattice for 
ML or the inversion symmetry between single layers for 
BL. These symmetries can be broken, either by adding 
hydrogen atoms to ML graphene^ ,2,31, or by a biased 
gate voltage applied to BL graphene [4]. The symme- 
try breaking is accompanied by opening of a gap in the 
spectrum. Then the question is whether or not such a 
gap is suppressed or supported by the Coulomb interac- 
tion. Previous studies have shown that disorder induces 
random fluctuations of the gap which can suppress the 
effective gap and allow ML and BL graphene to be a 
conductor and to have a metal-insulator transition for a 
sufficiently large average gap HQ. On the other hand, 
it has been discussed that a short-range (Gross-Neveu) 
electron-electron interaction can dominate the long-range 
Coulomb interaction, leading eventually to an insulating 
behavior Q. In contrast to this works, we will follow sub- 
sequently a more direct route to an insulator by assuming 
a small gap and study how this is affected by the Coulomb 
interaction itself. The problem of Coulomb interaction 
in graphene has been previously studied by employing a 
perturbative renormalization-group (RG) approach, for 
clean ML graphene [Ml as well as for disordered ML 
graphene |14| - [l6| . These studies show clearly a strong 
renormalization of the Fermi velocity in the clean case 
and considerable interplay between Coulomb interaction 
and disorder. 

The outline of this paper is as follows. In Sec. [Tl] we 
define the effective field theory for both graphene config- 
urations with Coulomb interaction. We introduce the 



gap into the action by hand and perform decoupling 
in the interaction channel by means of the Hubbard- 
Stratonovich transformation. We obtain expressions for 
the bare bosonic and fermionic propagators and interac- 
tion vertices. In Sec. Illll we write down renormalization- 
group flow equations for the gap parameter and fermionic 
wave- function renormalization factors and solve them for 
both gapless and gapped regimes. Furthermore, we ana- 
lyze the fixed points for both graphene configurations. 



II. THE MODEL 

We start with the zero temperature model for gapped 
ML and BL graphene. In the real space the Euclidean 
action of noninteracting ML and BL graphene in vicinity 
of a nodal point is given by 



^x- (1) 



X 



Here, Grassman fields V'^ — {i^AK'^BK'^BK'''^AK') 
represent four-component spinors on the sublattices A 
and B in the vicinity of nodal points K and K' in the mo- 
mentum space which depend upon the 2-1-1 dimensional 
vector X that contains imaginary time t and spatial vec- 
tor X as components. Matrices 3 = 1 (g) tr,- z = 1, 2, 



where (Jj 3 denote usual Pauli matrices. For ML graphene 
the operator V reads 



V = hvd^ 



(2) 



where v — \/3ta/2h denotes the bare (nonrenormalized) 
Fermi velocity and 9- usual differential operators. For 
BL graphene it has the components: 



fi2 



a: 



2/ij 



(3) 
(4) 



2 



with the bare band mass of electrons defined as 



3^2 a2 



Here, t and are in- and out-of-plane hopping energies 
respectively; a denotes the lattice spacing. The spectral 
gap Aq is simply introduced by hand. The instantaneous 
interaction is the same for both graphene configurations: 



hg 



X X' 



The microscopic strength of the Coulomb interaction be- 
tween electrons is given by 



ac 
Ye' 



where e denotes the elementary charge, ep the dielectric 
constant of the vacuum, a the fine structure constant, c 
the speed of light in vacuum and e the relative dielectric 
constant of the substrate. After performing a Fourier 
transform we obtain for both configurations 



4 



2 



ihqQ + h - T + AgTg 



with different kinetic energy parts. The integrals over 
momentum and frequency Q — (qq , cf) with the absolute 
value of the momentum q and zero-temperature Matsub- 
ara frequency Qq read Jq — {2it)^^ j dqQ(Pq and should 
be thought of being regularized by means of an UV- 
cutoff Aq. Furthermore we have re-scaled the interac- 
tion strength by the factor 27r that appears after Fourier 
transform, introducing g' = 2'Kg. The fermionic densities 
are defined as 



PQ 



For ML graphene the components of the vector h in 
the non-interacting part of the action read 



hvqi, 



while for BL graphene 

ft2 



hi = ^{qi -92), ^2 



-<ll<l2- 



(7a) 



(7b) 



Below we will assume h = 1. 

Now we map the pure fermionic action Eq. ([6]) onto the 
action containing both fermionic and bosonic degrees of 
freedom by means of the Hubbard - Stratonovich trans- 
formation [ITJ, 

S, [V-^ , V-] + ^] ^ 5o [V-t ,ij]+S,[c^]+Sy[i>\i:, cf,] ,(8) 



where the free bosonic action reads 

1 



2ff' 



(9) 



Q 



and the third term denotes the interacting Yukawa term 
describing coupling between fermions and bosons: 



(10) 



K 



From the mixed action Eq. we obtain vertices and 
propagators by taking functional derivatives with respect 
to each field. From 



S'S 



-{27r)%n>Go\Q), (11) 



we obtain the inverse fermionic propagator 

G^\Q)^iqa + h-f+Af,T3, (12) 



with components of the vector h defined in Eqs. (|7ap and 
(ITR. and from 



-{2^f5n,^Q'F-\Q) (13) 



(6) the inverse bosonic propagator 



Finally, the bare Yukawa vertex is obtained as 

8^S 



We arrive at 



r(Pi;P2,^3) = *(2^)''5p„p,+P3 



(14) 



(15) 



(16) 



III. RENORMALIZATION GROUP EQUATIONS 

The functional RG is conveniently defined in terms of 
the field dependent functional of effective action £[$], 
which in turn represents the Legendre transform of the 
generating functional of connected Green functions. For 
our purposes the ensemble average field $ = {'tp,^\(f>) 
is supposed to contain both fermionic and bosonic en- 
tries [13|- The functional £ depends on the IR-cutoff 
A ^ Aq, which is eventually removed. The derivation 
of the functional RG flow equation is described in detail 
in Refs. 17-20]. The RG flow of £ is generated by the 
regulator function introduced into the propagator of the 
non-interacting system and is determined by 



5A^A[1']--^Tr|aA[G„-],j(^ 



PCI 



, (17) 



3 



where [Gq ^ ] is the propagator of the non-interacting 
system depending on the cutoff A only via the regulator 
function. The matrix 



-{[Go 



-1 



(18) 



$=0 



denotes the regularized full inverse propagator with 
meaning the irreducible self-energy. The choice of the 
regulator will be specified a few lines below. Note that all 
quantities which appear on the right hand-side of Eq. (jl7p 
dwell on the space of composite fields $ and therefore 
represent 3x3 matrices. 

Since our main interest is the determination of 
the spectrum renormalization of fermions due to the 
Coulomb interaction, we will focus on the coupling pa- 
rameters in the fermionic sector of the theory. In the 
simplest approximation we make the following ansatz for 
the running effective action 



-2 / F'\Q)<j>Q4>_Q 



i^Qi^Q+K4>- 



(19) 



Q K 



which takes only the renormalization of the energy gap 
and of the electronic dispersion into account. We do not 
consider the renormalization of the Matsubara frequency 
since we assume the Coulomb interaction to be absolutely 
instantaneous. The inverse bosonic propagator F^^{Q) 
is defined in Eq. (ITil) . For momenta larger than the UV- 
cutoff Aq the action in Eq. ([T9| must reproduce the bare 
action from Eq. 



Therefore the initial conditions are 



chosen as =1 and A. 



Taking functional derivatives with respect to both 
Grassmanian fields on both sides of Eq. (fTT]) and putting 
subsequently $ = we arrive at the RG flow equation 
for the inverse renormalized fermionic propagator. For 
details of its derivation we reffer to Refs. |17lll8| . If we 
employ the regularization scheme with the regulator built 
in the fermionic lines only, this equation can be written 
in the following algebraic form (note an additional minus 
sign due to Fermi statistics): 



d^Gl\Q)= / G^{P)F{P-Q), 



(20) 



where F{Q) is the bare Coulomb potential defined in 
Eq. , and the single scale propagator is defined as 

Gyv = ^AlGp^ij^lGA (21) 

We will work within the so-called sharp cutoff regular- 
ization scheme |20] . Then the momentum cutoff is intro- 
duced as follows: 



Go.fl,(Q) = e(A<<z<Ao)Go(Q), 



(22) 



where e(A < g < Ag) = e(Ao-g)-e(A-g) ^ Qiq-h) 
as Aq — > oo. For momenta smaller than the UV-cutoff 
Aq, the flowing fermionic propagator G^{Q) is 



^1 



El[q) 



and hence the single-scale propagator j20] 

Mo - Z^^h ■ T - A^Ts 



Ga(Q) = <5(<z-A)- 



El{q) 



(23) 



(24) 



where we have introduced Ej^(q) 



e\{q) with 



the renormalized spectra of free fermions €j^{q) — ^vq 
for ML and e^C?) = (2Ai^A)"^9^ for BL. For both ML 
and BL the flow equations for the coupling parameter 
Ajy is extracted from Eq. (PU]) in the same way: 

dj,A^ ^ \Tr^{r,d^G^\Q)}\Q^^, (25a) 

where Trj denotes a trace operator acting on the pseu- 
dospin and valley space only. The RG flow equations 
for the factor are extracted differently for ML and 
BL due to the different scaling of the spectra in these 
configurations: 

ML : dj,Z^' = i|-Tr2 {T,d^Gl\Q)}\Q^^,{25h) 



BL 



5a^I' = J^Tr, {r.d^Gl^Q)}]^^^. 



(25c) 



for i = 1, 3^. Introducing the logarithmic flow parameter 
£ = log(AQ/A) we obtain the same flow equation for the 
gap for both graphene configurations 



d,A, 



gA^A 



71' 



(26a) 



where g — g' /Att^ but different equations for the wave- 
function renormalization factor: 



ML 



BL 



1 -gZA 



3 gZ^K 



71' 



(26b) 

(26c) 



Here we have used the identity d^^Z^^ = —Z'^'^d^Zg. The 
flowing free fermion spectra are 

= vZ'[^K, 



ML 
BL 



(2mZ, 



(27a) 
(27b) 



The scaling dimension of the energy, (i.e. the dynamical 
exponent) is defined as z = 1 — ry^ for ML and z = 2 — 
for BL. Here, rj^ is referred to as the anomalous dimension 
which can be obtained from the parameter Z^ by 



-dglogZg 



(28) 



Below we discuss solutions of these equations in gapless 
and gapped regimes. 
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A. Gapless regime 



In the gapless regime Eqs. (|26a|) - (|26c|) are easily solved. 
The only solution of Eq. (|26ap is the trivial one = 
^£=0 — 0' while Eqs. (j26bp and p6cp reduce to 



A 



ML 



with Aj^^L — d/i'^'^) foi' and correspondingly 
with AgL = 3/i.g/(4AQ) for BL with solutions: 



ML 
BL 



— 1 + A^L^, 



1 - ^BL + ^BLe^- 



(29) 
(30) 



The result of Eq. ([29)) corresponds to the well- 
known logarithmic renormalization of the Fermi velocity 
= Z^^v in clean ML graphene due to the Coulomb 
interaction [i-[ll|l2l|. Similarly, Eq. §UI describes the 
renormalization of the electronic band mass /x^ = Z^ii in 
BL. At small momenta the band mass decreases propor- 
tionally to the momentum /i^ cx A, i.e. the particles be- 
come effectively faster in analogy to ML. Using Eq. 
we obtain expressions for the anomalous dimension 



ML : 77^ = Xml^i, 
BL : 7jf = AgL^^e^. 



(31a) 
(31b) 



For £ ^ Eq. (|31ap approaches zero, meaning that the 
scaling dimension of the energy in ML remains z = 1 and 
nothing changes the relativistic behavior of electrons. In 
contrast, Eq. pibl) approaches in this limit unity. This 
means that the scaling dimension of the energy z — 2 — 
becomes relativistic in BL with the velocity — 85/8, 
i.e. in vacuum c/v^ w 1450. Therefore, in absence of 
a gap in the spectrum of BL the Coulomb interaction 
attempts to linearize the fermionic dispersion in vicinity 
of the nodal points. Similar conclusions have been re- 
cently made by Kusminskiy et al. (23)] for finite values of 
chemical potential. Their findings provided a good ex- 
planation of recent cyclotron experiments [24| , where the 
effects discussed here have been observed. 

An estimation for the suitable scale below which this 
effect is observable can be made as follows: The only 
scale which affects the flow of the band mass in gapless 
BL graphene can be read off from Eq. ([50]) (cf. Fig. S]): 



log 



1 - A 



BL 



A 



BL 



(32) 



Choosing Aq to be equal to the inverse lattice spacing 
we find for a realistic experimental situation (e = 1 4) 
£' « 2.3 3.8 and the corresponding momentum scale to 

-2^7- 10-2 a"\ 



be of the order 



An 



1 • 10- 



A, /An 



FIG. 1: Renormalization of the gap in ML because of 
Coulomb interaction. The crossover scale is determined 
from Eq. ((35}. The initial value of the gap is Aq = 0.2woAo, 
the dielectric constant e = 1. The dashed line shows the 
Kane/Mele asymptote from Eq. (|34p . The crossover scale I, 
is determined from Eq. p5|) . 
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FIG. 2: Renormalization of the Fermi velocity vi = VqZ'[^ 
in ML because of Coulomb interaction with (solid line, full 
solutions of Eqs. (I26ap and (|26b|l l and without a gap (dashed 
line, Eq. ((29])). The crossover scale 1^, is determined from 
Eq. ((35)). 



B. Gapped regime 

Naively, for ^ we can neglect A^ in the denom- 
inator. For ML we obtain from Eqs. ()26ap 



(33) 



which together with Eq. (|29p reproduces the Kane/Mele 
result pj,]: 

A, = Ao(l-f Aml^)' = AoV, (34) 

with the apparently logarithmically diverging gap. How- 
ever, Eq. (|34| suggests that at large i denominators in 
Eqs. (|26ap and (|26bp are dominated by the gap, i.e. 

? ~ A«, where the crossover scale 



E, 



can be determined from the condition 



A, 



(35) 



which turns out to be a nonlinear algebraic equation if 
we take Eqs. (|27ap . (|29|) and (|33|) into account. However, 
Eq. psp can be uniquely solved numerically. The solution 
of Eq. (|26ap in this case becomes 



A, « A, +5A,(1 



(36) 





FIG. 3: Renormalization of the gap in BL because of 
Coulomb interaction. The crossover scale is determined 
from Eq. ((35}. The initial value of the gap is Aq — 0.2woAo, 
the dielectric constant e — 1. The dashed line shows the large 
kinetic energy asymptote from Eq. (|38|l . The crossover scale 
l^, is determined from Eq. H35[). 



FIG. 4: The renormalization of the band mass fie = fio^e due 
to the Coulomb interaction. The solid line shows the flow of 
the band mass of BL graphene with a gap. Dashed line shows 
asymptotic renormalization without a gap. In this case the 
electronic band mass scales to zero. This leads to the linear 
scaling of the spectrum. The scale £' is found from Eq. (f32)l. 

The physical gap is obtained for £ ^ oo 

A, = A, +.gA, +gA,. (37) 

Therefore the Coulomb interaction in ML supports the 
gap once it is opened, independently of the bare gap mag- 
nitude. A typical flow of the gap parameter in ML is 
shown in Fig.[TJ At the same scale the logarithmic growth 
of the Fermi velocity stops and it also stabilizes at the 
finite value « + gA^/A, as depicted in Fig. [2] 

The solutions for the gap in BL are similar in spirit 
but with an extra fixed point. For A^ ^ we obtain 

A, = AoZ7«/3 (38) 

and therefore Ai^^o oo. The solution for A^ ^ 
is formally given by Eq. (l36t . too, such that the fiow of 
the gap stabilizes at some finite value A, (cf. Fig. [3]). 
On the other hand, the presence of the gap stabilizes the 
flow of the wave function renormalization factor and 
therefore the flow of the electronic band mass /i^ — (J-qZ^, 
which in this case remains finite (cf. FigH}. The scaling 
of the kinetic energy is in this case also preserved and 
remains equal to 2. 



FIG. 5: Schematic RG flow for both graphene configurations 
in the space spanned by the dimensionless kinetic energy 
and gap parameter A^. 

In order to shed some light on the topological prop- 
erties of the RG flow in the parametric space it is con- 
venient to redefine Eqs. (|26bp and (|26cp in terms of ki- 
netic energy and introduce dimensionless parameters by 
expressing both the gap and kinetic energy in units of 
Coulomb energy: 

A, . |, (39a) 
= (39b) 

with defined in Eq. ([27a| for ML and in Eq. (f27al) for 
BL. For both ML and BL we arrive at the same equation 
for the rescaled gap 

5,A, = A,-f ^^4=, (40a) 

while equations for the rescaled kinetic energy are differ- 
ent due to different scaling behavior of spectra: 

ML: d,e,^l^J^, (40b) 

BL: d,e, = l^jL=--e„ (40c) 

The flow in the parametric space is schematically shown 
in Fig. ([5]). The fixed points (FPs) are obtained by set- 
ting the right-hand sides of Eqs. (|40ap - ()40cp to zero and 
solving the emerging system of algebraic equations. For 
ML graphene the only instable fixed point is at both 
A^ = and = 0. From Eq. (pQb} follows that this 
fixed point can be reached if 

= ^ ^ 0. (41) 
9 

Since Z^^ flows to a flnite value, this can only be sat- 
isfled if 5 — >■ cxD. This is a case of the famous quantum 
phase transition discussed in Refs. (8HTl|. The instabil- 
ity of the fixed point means that the flow can leave it 
in every direction. For any finite initial value of the gap 



6 



it develops infinitely large value which indicates a finite 
physical gap. In contrast to the ML graphene, there is a 
nontrivial fixed point at finite dimensionless kinetic en- 
ergy = 3/8 in the case of gapless BL graphene. This 
fixed point is characterized by the anomalous scaling di- 
mension rj^ — 1, i.e. the spectrum of BL becomes in this 
case linear. However this fixed point is instable with re- 
spect to the finite gap direction, i.e. once a small gap is 
opened in the spectrum the flow cannot reach this fixed 
point anymore but runs towards an infinite value. On the 
other hand, since the numerical value of at this fixed 
point suggests a strong coupling regime we might need to 
go beyond the truncation Eq. ([T9)) and take additionally 
flow of the tjjtjj'^ (j)-vej:tex into account. 



IV. CONCLUSIONS 

In conclusion, we have studied both ML and BL 
graphene with Coulomb interaction and a uniform gap 
by employing a renormalization-group technique. In con- 
trast to previous approaches to gapped ML graphene 
based on the renormalization group approach [? lol flllfT^ . 
which predict logarithmically divergent renormalization 
of the gap and the Fermi velocity, our results suggest a 
saturation of RG flows at an intrinsic scale related to the 
gap. This saturation takes place for both ML and BL 
graphene, for any finite initial value of the gap no matter 
how small it is, and since measured quantities should be 



flnite, this might be suggestive of a gap in the spectrum 
of both configurations at energies below 0.1 eV. 

In ML graphene the Coulomb energy exhibits the same 
scaling as the kinetic energy. Once a spectral gap is 
opened it creates an additional length scale which dom- 
inates the physics at small momenta. This scale cuts off 
the logarithmic divergence of the Fermi velocity and gap 
itself such that the flow of both quantities stabilizes at 
the finite value. 

For gapless BL graphene is shown that Coulomb inter- 
action renormalizes the electronic band mass which scales 
to zero for small momenta. This leads to a paradoxical 
result that the electronic spectrum should become linear 
close to the charge neutrality point. This regime cor- 
responds to a stable fixed point and therefore the flow 
should inevitably go into this point. The quadratic scal- 
ing of the spectrum is rescued by the presence of the gap, 
since for any finite starting values of the gap the flow of 
the band mass always saturates at a finite value. 
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